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We present a new model for the formation of spherically symmetric clusters in an ex- 
panding Universe. Both the Universe and the collapsing cluster are governed by the same 
pressure less fluid equations for which a uniform initial density profile is assumed. A sim- 
ple perturbation imposed on the initial velocity field gives rise to an over-density which 
closely models real clusters. The computation of photon paths allows us to evaluate the 
gravitational effects imprinted on a Cosmic Microwave Background (CMB) photon passing 
through such an evolving mass. We also consider the lensing properties of collapsing clus- 
ters and investigate the effect of a population of such clusters on the primordial microwave 
background power spectrum. 
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1 Introduction 

There exist several mechanisms which produces anisotropics on the CMB: the primary anisotropics 
which take place at the last scattering epoch and the secondary anisotropics such as the Sunyaev- 
Zel'dovich effect which occurs while the CMB photons are travelling through the universe. We will 
concentrate here on a particular secondary anisotropy first discussed by Rees & Sciama (1968). As 
pointed out by these authors, there is a gravitational effect on the CMB photons while they are cross- 
ing evolving cosmic structures such as a collapsing cluster of galaxies or expanding voids. Indeed 
photons climb out of a slightly different potential well than the one into which they entered. This 
effect is usually referred as the Rees-Sciama effect. 

Early approaches to this problem were based on the "Swiss Cheese" model (Rees & Sciama (1968), 
Dyer (1976), Kaiser (1982), Nottale (1982) & Nottale (1984)), whereas more recent attempts have 
used the continuous Tolman-Bondi solution (Panek (1992), Saez et al. (1993)). Our aim is to present 
here an improved model which, for the first time, treats such a problem exactly (see section 

In section || we apply our model to a rich galaxy cluster at a redshift z = 0.08. We investigate 
the properties of the collapsing cluster and verify that its density-profile is of a realistic shape. We 
then give, in section ^, an upper limit for the temperature decrement that such an evolving cosmic 
structure could imprint on the CMB. Our result is compared to previous works. 

Section 1^ deals with distant clusters (z = 1) since it is reasonable to assume that, at such an epoch, 
those objects were in a stage of formation and that our spherical collapse model would be suitable. 
The focus is put on the gravitational lensing phenomenon and how this would distort the observed 
CMB power spectrum. 

We finish in section || by a short conclusion and discuss the applicability of our work to model the 
microwave decrement observed towards the, possibly lensed, quasar pair PC1643+4631 A&B. 



The methods and results introduced in this poster are discussed in further details in forthcoming 
papers (Lasenby et al. (1997), Dabrowski et al. (1997)). 
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2 Model 



In this poster we are concerned with the formation of spherical galaxy clusters (i.e. pure infall) in an 
expanding background universe. As in the Tolman-Bondi model applied by Panek (1992), the evolu- 
tion of both the collapsing cluster and the expanding Universe are governed by the same pressure-less 
fluid equations. Initially, the density of our fluid is uniform. Only the velocity field is perturbed in 
order to form a realistic cluster. In a given cosmology, this perturbation is controlled by two parame- 
ters (ro and a) as illustrated in figure |^. tq is related to the size of the perturbation and a to its rate 
of growth. Those parameters are free and we choose them so that the resulting observed cluster is 
realistic (see section ||) . 




radius 



Figure 1: Fluid initial conditions for our model. Two parameters control the perturbation: ro and a. 



We find that our model has many advantages and treats the problem much more rigorously than 
older attempts, particularly in the case of works based on the "Swiss Cheese" (SC) models. Indeed, 
as shown in figure ^, the SC models deals with discontinuous density and velocity distributions since 
an unrealistic vacuum region is needed to separate the perturbation from the Universe. Another 
improvement made with regard to the SC models is that our perturbation is defined by only two free 
parameters rather than three in the SC case. 




Figure 2: Fluid initial conditions for Swiss-Cheese models. Three parameters control the perturbation: 
Pc ri and a. 

More recent studies of the problem use the Tolman-Bondi solution (Panek (1992), Saez et al. 
(1993)), which gives a similar approach to our own. In particular, both the forming cluster and 
the external Universe are treated as a whole and the fluid distributions are continuous. We believe 
our model is a refinement of the work based on the Tolman Bondi solution. Our mathematical 
formulation is clearer and does not need any approximation. The underling gauge theory of gravity 
(Lasenby, Doran & Gull (1997)) that we use deals straightforwardly with observable quantities such 
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as the photon energy (Lasenby et aL (1997)). Therefore, we do not need to remove a posteriori dipole 
nor quadrupole-like anisotropy form the obtained CMB fluctuation (e.g. Panek (1992)). 

3 Cluster properties 

In this section, we model the formation of a very rich Abell cluster. The characteristics (see table ^ 
are chosen to be similar to those of the clusters described in Quilis et al. (1995) and Panek (1992) so 
that our results can be compared. 



Table 1: Cluster characteristics. 





Ho ^50 


Ho = 100 


Redshift 


0.08 


0.08 


Distance 


450 Mpc 


225 Mpc 


Max density 


0.7x 10*pm-^ 


1 X lO^pm^"^ 


Core radius 


0.46 Mpc 


0.23 Mpc 


Mass r < 4 Mpc 


1.3x lO^^'M© 


5.8 X IO^'^Mq 



The cluster is placed at a redshift z — 0.08, with a core radius Rc = 0.23 /iq ^ Mpc (i?c is defined as 
the radius at which the cluster energy falls to one- half its maximum value). The maximal baryonic 
density Pmax encounter by an observed photon is taken to be lO'* /iq^^ protons m~'^ and the baryon 
to dark-matter mass ratio is assumed to be 0.1. The /i-dependencies for Rc and Pmax ensure that 
the observed cluster characteristics (i.e. angular size and X-ray luminosity) are independent of the 
Hubble parameter. Throughout this paper we assume Qq — 1. 

Figure ^ shows the fluid density and the fluid velocity as a function of proper time as experienced 
by a photon which travels straight through the centre of the cluster. From the figure we see clearly 
that the fluid distributions are continuous. 




Figure 3: Plot of the fluid density (left) and fluid velocity (right) as experienced by the photon while 
it is travelling through the Universe and the collapsing cluster. The horizontal axis marks the time 
prior to the present epoch in miUions of years. We assume ho — 0.5. 

We next consider (fig ^ the density profile of the cluster at the time tc = 1436.1 Myr ago, when the 
photon experienced the maximum density 0.7 x 10^ protons m^"^. In order to check that the obtained 
density distribution is of realistic shape, we fit our profile to the equilibrium spherical King model 
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(King (1966)): 



P[r) = Pmax 




where pmax is the central density of the cluster, R represents some core radius and /3 is the conventional 
power-law index for such models. The results of a simple least-squares fit are pmax — 0.73 ±0.05 x 
10'' protons m"'^, i? = 0.52±0.05 Mpc and /3 = 0.8±0.15. The result of this fit is very encouraging since 
it shows that our model leads to a radial density profile that matches quite closely those of observed 
clusters. 




Figure 4: Plot of the density when the photon is at the centre of the cluster. The obtained profile is 
realistic since it can be fitted by a spherical King model. 



4 Effect on the CMB temperature 

We now consider photon energy and concentrate on the CMB anisotropy produced when the photon 
passes through the cluster described in section ^. One must remember however that our model as- 
sumes a pressure-less fluid and so the cluster, in pure-infall, might evolves too rapidly. Therefore the 
temperature perturbations given in this section should be considered as upper limits. 

Figure || shows the CMB anisotropy due to the gravitational perturbation of the cluster described 
m table |l|. The maximum temperature distortion occurs at the centre of the cluster and has the value 
AT/T = -0.96 X 10"^ and AT/T = -5.2 x 10"^ for /iq = 1 and /iq = 0.5 respectively. One notices that 
the distortion extends to rather large projected angles (e.g. AT/T ^ — 1 x 10^^ for an observed angle 
of ~ 2.5 degrees in the /lo = 0.5 case). We may also point out the fact that the anisotropy becomes 
slightly positive at large angles. 

We can compare the central decrements calculated above with those of previous authors. For 
Panek (1992) type I and type II clusters, the calculated central decrements are AT/T = —1.5 x 10^^ 
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Figure 5: Temperature perturbation l^T/T imprinted on the CMB as a function of the observed 
angle from the centre of the cluster. 



and AT/T — —6.0 x 10"^ respectively, whereas Quilis et al. (1995) quote the decrement /\T/T — 
— 1.2 X 10^^. For two cluster models with similar physical properties Chodorowski (1991) finds /^T/T = 
—7.7 X 10^^. Finally Nottale (1984) used the SC model to predict a considerably larger central 
decrement of AT/T ^ 10^'*. However, this last value corresponds to a very dense, unrealistic cluster. 
We can notice that, for ho = l, our predicted value is in rough agreement with Panek (1992) type II, 
Chodorowski (1991) and Quihs et al. (1995). However our result for Hq = 0.5 is about five times larger 
since previous works predict the same value whatever is ho- We then suggest that, in such a Universe, 
the effect on CMB photons may be more significant that previously stated. 

5 Distant clusters 

We consider, in this section, the formation of a massive cluster located at a redshift z — 1. It is 
more reasonable to apply our model for a distant cluster rather than the one described in section |[ 
Indeed, at such redshifts one could expect galaxy clusters to be in a formation process and display large 
inward radial velocities. As from now, the cosmology is taken to be fio = 1 and Hq ~ 50kms~^Mpc~^. 
The cluster characteristics are the following: z = l, maximum density = 1.0 x 10'*pm~'^, coreradius = 
0.46 Mpc and the mass within 4Mpc is 1.9 x 10^^ Mq. In the next section we investigate the lensing 
effect caused by such a cluster as well as the effect imprinted on the CMB temperature and power 
spectrum. 

5.1 Gravitational lensing and effect on the CMB 

A massive cluster such as described above is a powerful gravitational lens. This effect is illustrated in 
figure ^. Our model also allows a quantitative study of the effect of dynamically evolving lenses on the 
CMB fluctuations. By simulating maps of fluctuations in the CMB due to inflation and/or topological 
defects, we can investigate how features in these maps are affected by the presence of the lens. As an 
example, flgure shows the effect of our collapsing cluster on a CMB fluctuations realisation. The 
result is a central decrease of temperature at the centre (i.e. 180 /iK in our case) as well as a global 



5 



magnification of tlie features. 




Figure 6: Lensing image of 4 sources at a redshift z — 3.8 by the rich cluster described in section ^ 
The blue bullets (A,B,C,D) show the positions of the sources as seen in the absence of the cluster. 
The blue cross marks the centre of the cluster. The source A placed on the lens caustic displays two 
images, one radially and the other tangentially elongated. 

Aside from studying individual maps, we can also study the statistical effect that a population 
of such lenses would have on the power spectrum of the CMB fluctuations. This power spectrum 
is currently of immense theoretical and observational interest, since it is now a possibility that the 
spectrum may be measured, and so the values of cosmological constants may be found to unprecedented 
accuracy. Such determinations depend, however, on distinctive features in the power spectrum which 
may be affected by a population of large clusters. Previous work on this effect has been performed 
such as Seljak (1996), but we have now incorporated the effect of dynamically collapsing clusters in 
a proper general relativistic fashion. Figure || gives the result obtained with the cluster used in the 
current section. 

We find that the distinctive peaks in the power spectrum are slightly smoothed out by such a 
population, and so this effect should be taken into account in future analyses. 

6 Discussion 

A CMB decrement of ~ 380 /i Jy has been observed by the Ryle- Telescope towards the pair of quasars 
PC1643+4631 A & B at red-shift ^ 3.8 separated by 100 arc-seconds (Jones et al. (1997)). Assuming 
a Sunyaev-Zel'dovich effect from a cluster of galaxies, this is indicative of a rich intervening cluster 
of total mass ~ 10"'^^ M0, although, no X-ray cluster has yet been observed in that direction. This 
suggests that the object is lying at a red-shift greater than one. The quasars' spectra analysis reveals 
that they might be originating from a unique source lensed by the distant cluster (Saunders et al. 
(1997)). 
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View without the cluster 



with the cluster 





Figure 7: Effect of the collapsing cluster described in section ^ on the CMB fluctuations. The first 
map is a simulation of the CMB fluctuations due to inflation. The second one represents the same 
patch of the sky (4x4 degrees) but with a rich cluster of galaxies at its centre. 









Initial PS 

• Lensed PS 




Figure 8: Effect on the CMB power spectrum of the rich cluster described in section |[ We assume 
one cluster per 10x10 degrees field. The red curve is the unperturbed power spectrum as the black 
dots represent the "lensed" spectrum. 
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Using our model, we are able to fit those observations and retrieve the S-Z flux decrement as well 
as the 100 arc-seconds separation. To do so, we have to consider the cluster described in section || 
with an electron temperature of 2 x 10^ K, which corresponds to a rich but rather cool cluster. As 
we stated in section the resulting Rees-Sciama temperature decrement is < 180 which is a 
considerable fraction of the S-Z decrement, here ~450/iK. Fitting those data is of great importance 
since such a distant lensed-pair has previously never been observed. As the cluster lies at a rather 
large red-shift, we believe that a model such as ours, computing rigorously both the cluster and the 
universe evolutions, is needed. 



References 

Chodorowski M., 1991, MNRAS, 251, 248. 

Dabrowski Y., Hobson M.P., Lasenby A.N., Doran C.,1997, MNRAS, submitted. 
Dyer C.C., 1976, MNRAS, 175, 429. 

Jones M.E., Saunders R., Baker J.C., Cotter G., Edge A., Grainge K., Haynes T., Lasenby A., Pooley 

G., Rottgering H., 1997, ApJ, 479, LI. 
Kaiser N., 1982, MNRAS, 198, 1033. 
King LR., 1966, Astron. J., 71, 64. 

Lasenby A.N., Doran C.J.L, GuU S.F, 1997, To appear in: Phil. Trans. R. Soc. Lond. A. 
Lasenby A.N., Doran C., Hobson M.P., Dabrowski Y., Gull S.F, 1997, MNRAS, submitted. 
Nottale L., 1982, A&A, 110, 9. 
Nottale L., 1984, MNRAS, 206, 713. 
Panek M., 1992, ApJ, 388, 225. 

Quilis v., Ibahez J.M., Saez D., 1995, MNRAS, 277, 445. 

Rees M.J. & Sciama D.W., 1968, Nat, 217, 511. 

Saez D., Arnau J.V., FuUana M.J., 1993, MNRAS, 263, 681. 

Saunders R., Baker J.C., Bremer M.N., Bunker A.J., Cotter G., Eales S., Grainge K., Haynes T., 

Jones M.E., Lacy M., Pooley G., Rawlings S., 1997, ApJ, 479, L5. 
Seljak U., 1996, ApJ, 463,1. 



8 



